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A promising approach to solving hard binary 
optimisation problems is quantum adiabatic an¬ 
nealing (QA) in a transverse magnetic field. An 
instantaneous ground state — initially a symmet¬ 
ric superposition of all possible assignments of N 
qubits — is closely tracked as it becomes more 
and more localised near the global minimum of 
the classical energy. Regions where the energy 
gap to excited states is small (e.g. at the phase 
transition) are the algorithm’s bottlenecks. Here 
I show how for large problems the complexity 
becomes dominated by O(logA) bottlenecks in¬ 
side the spin glass phase, where the gap scales 
as a stretched exponential. For smaller A, only 
the gap at the critical point is relevant, where it 
scales polynomially, as long as the phase transi¬ 
tion is second order. This phenomenon is demon¬ 
strated rigorously for the two-pattern Gaussian 
Hopfield Model. Qualitative comparison with the 
Sherrington-Kirkpatrick Model leads to similar 
conclusions. 

Quantum algorithms offer hope for tackling com¬ 
puter science problems that are intractable for classi¬ 
cal computer^ However, exponential speed-ups seen 
in, e.g. number factoring, have not materialised for 
more difficult NP-complete problem^. Those problems 
are targeted W the quantum adiabatic annealing algo¬ 
rithm (QAplEl. Any NP-hard problem can be recast as 
quadratic binary optimisation. QA solves it by imple¬ 
menting a quantum Hamiltonian, written with the aid of 
Pauli matrices as 

N N N 

-^ = “2 X - X ~ X (1) 

i,k=l i=l i=l 

Here the first two terms, diagonal in z-basis, encode the 
objective function. The last term represents the mag¬ 
netic field in the transverse direction, which is decreased 
from r(0) ^ 1 to r(Tann) = 0- The time Tann needed 
by the algorithm is determined by a condition that the 
annealing rate is sufficiently low to inhibit non-adiabatic 
transitions: 

dr/dt<AA.Ar. ( 2 ) 

These are most likely near points where the instantaneous 
gap to excited states AE attains a minimum as a function 
of T; further, AT is defined as the width of the region 
where the gap remains comparable to its minimum value. 

QA offers no worst-case guarantees on time 
complexit}i3, but initial assessments of typical case 


complexity were optimistic. Both experimentaP and 
theoreticaP evidence hinted at performance improve¬ 
ment over simulated annealing for finite-dimensional 
glasses; however, some empirical evidence in support of 
the theory has recently been called into questiorP^. Early 
exact diagonalisation studies also observed polynomially 
small gaps in the constraint satisfaction problem on 
random hypergraph^^, but that finding had been chal¬ 
lenged by quantum Monte Carlo studies involving larger 
size^^. Benchmarking of a hardware implementation of 
QA, courtesy of D-Wave Systems, shows no improvement 
in the scaling of the performanc d^^ l ^^ i Whether that 
might be attributable to a finite temperature at which 
the device operates or its intrinsic noise is unclear at 
presenlP^l^. 

Statistical physics offer some intuitive guidance: Small 
gaps develop at the quantum phase transition point and 
become exponentially small when the transition is 1®^- 
ordeiC^HlIl 

or when different parts of the system become 
critical at different times for strong-disorder continuous 
phase transitions^. The most promising candidates for 
QA are thus problems with bona fide 2^^-order phase 
transitions, where the disorder is irrelevant at the QCP. 

The scaling analysis described here suggests a polyno¬ 
mially small gap at the critical point of the archetypal 
spin glass: the Sherrin gton- Kirkpatrick (SK) modeP^^^. 
It has been pointed ouiP^ that QA may still be doomed 
by the bottlenecks in the spin glass phase. Exponen¬ 
tially small gaps away from the critical point have been 
observed in simulations^, but adequate theoretical de¬ 
scription of this phenomenon has proven challenging. A 
perturbative argument in support of this qualitative pic¬ 
ture has been considered in ref. [561 However, the re¬ 
sults were not borne out by more accurate analysis that 
took into account the extreme value statistics of energy 
levelS^. 

The present manuscript sheds new light on the mecha¬ 
nism of tunneling bottlenecks in the spin glass phase. Us¬ 
ing exact, non-perturbative, methods, this is illustrated 
for a simple model, but the main findings are expected 
to be valid for quantum annealing of more realistic spin 
glasses. During annealing, the system must undergo a 
cascade of tunnelings at some specific values of Ti, r 2 ,... 
in an approximate geometric progression. Eor a finite sys¬ 
tem size, these bottlenecks are few, O(logA), and may 
not even appear until N is sufficiently large, highlight¬ 
ing the challenge of interpreting the results of numerical 
studies. Bottlenecks also become increasingly easier as 
T ^ 0, counter to expectations that tunnelings are in¬ 
hibited as the model becomes more classical. A related 
finding is that the time complexity of QA is exponential 
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only in some fractional power of problem size: a mild 
improvement over more pessimistic estimateJ^. 


RESULTS 

Summary. The spin glass phase, which is entered 
below some critical value of the transverse field Pc, is 
characterised by a large number of valleys. Often, this 
transition is abrupt, driven by competition between an 
extended state and a valley (localised state) with th e low- 
est energy, as occurs in the random energy modeP^^. 
The exponentially small overlap between the two states 
then determines the gap at the phase transition. How¬ 
ever, even if new valleys develop in a continuous manner 
as r decreases, small changes in the transverse field may 
result in a chaotic reordering of associated energy levels, 
leading to Landau-Zener avoided crossings and attendant 
exponentially small gaps. 

Nonetheless, attempts to make this intuition exact are 
fraught with potential pitfalls. For increasing F, two 
randomly chosen valleys are equally likely to come ei¬ 
ther closer together or further apart in energy. In the 
case of the former — and further, if the sensitivity of en¬ 
ergy levels to changes in the transverse field is so large 
that the levels ‘collide’ before either valley disappears — 
avoided crossing will occur. This may not be necessarily 
the case when one considers ‘collisions’ with the ground 
state, which are of particular concern to QA. The ground 
state corresponds to a valley with the lowest energy; this 
and other low-lying valleys obey fundamentally different 
statistics of the extremes. 

A case in point is the analysis of ref. [26j which devel¬ 
ops perturbation theory in F. The classical limit (F = 0) 
is used as a starting point; how that analysis might be 
extended to F > 0 has also been discussecP^. A type 
of constraint satisfaction problem (CSP) has been con¬ 
sidered: classical energy levels are discrete non-negative 
integers (number of violated constraints) but have expo¬ 
nential degeneracy. ‘Zeeman splitting’ for F > 0 scales as 
which, for large problems, may be sufficient to over¬ 
come the 0(1) classical gap and cause avoided crossings 
of levels associated with different classical energies. Yet 
this trend disappears if only levels with the smallest ener¬ 
gies (after splitting) are considered; these are relevant for 
avoided crossings with the ground state. This about-face 
is not immediately apparent, only coming into play for 
N > 100, when the exponential degeneracy of the classi¬ 
cal ground state sets in. It has, however, been confirmed 
with analytical argument and numeric^. Consequently, 
arguments based on perturbation theory cannot be used 
to establish the phenomenon. 

This manuscript offers a fresh perspective, illustrated 
by studying quantum annealing of the Hopfield model. 
Mean-field analysis correctly describes thermodynamics 
if the number of random ‘patterns’ is small. The method 
is further extended to extract information about exact 
quantum energy levels. Importantly, the classical energy 



FIG. 1. Sketch of the behavior of the gap in a Hopfield model 
with the Gaussian distribution of disorder variables: in the 
paramagnetic phase (F > Fc), in the glassy phase (F < Fc) 
and in the critical region (F Fc). Scaling of the typical 
gap in these regions is indicated in bold letters using big-0 
notation. The area F < Fmin ~ 1/N is where the discrete 
nature of the energy landscape becomes manifest: the ground 
state becomes nearly completely localised. The glassy phase 
also contains log A isolated bottlenecks, indicated with red 
arrows, where the gaps scales as a stretched exponential. 


landscape is made much more complex by insisting that 
the distribution of disorder is Gaussian. Fig. [^sketches a 
‘phase diagram’ obtained for this model. For decreasing 
F the gap changes as follows: (1) it is finite (does not 
scale with N) in the paramagnetic phase, F > Fc; (2) 
scales as in the narrow region of width 

around F = Fc; (3) increases slightly, with typical values 
scaling as for F < Fc. In addition, there are 

avoided crossings at isolated points Fi, F 2 ,..., which is a 
manifestation of the transverse field chaos, demonstrated 
conclusively in this work. 

The first non-trivial example requires two Gaussian 
patterns. In this case the ‘energy landscape’ is effec¬ 
tively one-dimensional, which greatly simplifies the anal¬ 
ysis. The most important element of this analysis is ac¬ 
counting for the extreme value statistics associated with 
the valleys (local minima) having the lowest energies. To 
this end, the distribution of energies of the classical land¬ 
scape must be conditioned so that they are never below 
the energy of the global minimum. This becomes feasi¬ 
ble when reformulated as a continuous random process, 
in the limit N ^ 00 . 

This particular model is not as interesting from the 
computer science perspective, not least because it affords 
an efficient classical algorithm. It is sufficiently simple so 
that a complete quantitative analysis presented later on 
in the manuscript has been possible. Yet, the model cap¬ 
tures the essential properties of the spin glass: its qualita¬ 
tive features directly apply to much more general models, 
including Sherrington-Kirkpatrick. The most important 
feature of the classical energy landscape is that exhibits 
fractal properties, which both ensures that hard bottle¬ 
necks are present in the spin glass phase and also governs 
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their distribution. The role of the transverse field is to 
approximately coarse-grain it on scales determined by T, 
eliminating small barriers; thus the number of valleys 
decreases as T grows. A specific random process, cor¬ 
responding to the energy landscape of the dnfinite’-size 
instance, will contain every possible realisation of itself at 
some ‘length scale’. Some realisations will contain high 
barriers that cannot be easily overcome; these will lead 
to tunneling bottlenecks. 

This intuition can be used to immediately establish 
the scaling of the number of tunneling bottlenecks. Since 
the model contains no inherent length scale in the limit 
A" ^ oo, it can be argued that the expected num¬ 
ber of tunneling bottlenecks in a finite interval [ri;r 2 ] 
should be a function of the ratio r 2 /ri. The loga¬ 
rithm is the only function that respects additivity, i.e. 
Munn([ri,r 2 ]) =Munn([ri;r])+Munn([r',r 2 ]). To ob¬ 
tain the total number of bottlenecks, one considers the 
interval [rniin;rc]. Here Tc ^ 1 is the boundary of the 
spin glass phase. The lower cutoff, Tmin, corresponds 
to the lowest energy scale of the classical model, which 
scales as an inverse power of system size, e.g. as 1/N 
for the Gaussian Hopfield model. In a sense, tunneling 
bottlenecks are connected to the T = 0 ‘fixed point’ (note 
that the classical gap vanishes asymptotically). To sum¬ 
marize, the number of tunneling bottlenecks will grow 
as 


Munn^alnA^. (3) 

Locations (in F) will depend on specific disorder realisa- 
tion, but self-similarity ensures that the successive ratios 
Tn/Tn-\-i converge to a universal distribution. 

This logarithmic rise is far weaker than a power law 
seen in some phenomenological models of temperature 
chao^ and, as has been argued above, likely to be a 
universal feature. The prefactor is model-dependent; its 
numerical value can be used to estimate the minimum 
problem size for which the mechanism becomes relevant, 
via A/tunn ~ 1- A value of a ^ 0.15 is obtained for the 
problem at hand, so that additional bottlenecks become 
an issue for N > 1000. Prior numerical studies similarly 
required large sizes before the exponentially small mini¬ 
mum gap was observecP^, and so far there has been no 
evidence of two or more exponentially small gaps coex¬ 
isting. The slow, logarithmic increase of the expected 
number of bottlenecks is the most likely culprit. 

A notable feature of these results is that tunnelings 
become progressively ‘easier’ as T 0 despite the fact 
that the model becomes more classical. Tunneling gaps 
increase as 


AA 


(n) 




( 4 ) 


Notice that they cease to be exponentially small for T < 
Tmin; at that point the ground state is already localised 
near the correct global minimum. The power law expo¬ 
nent for this stretched exponential is model-dependent, 
related to the scaling of barrier heights. These scale as 


N^/‘^ for the Gaussian Hopfield model, which, together 
with 0{N) scaling of the effective mass, gives rise to the 
^3/4 iqyir in the exponent. 

The finding that the gaps increase for smaller T de¬ 
serves explanation. Typically, valleys with similar ener¬ 
gies differ by up to N/2 spin flips. This changes once 
lowest energies are considered: All spin configurations 
with energies less than e above the global minimum are 
contained in a neighborhood of radius 0[(Ae)^/^], us¬ 
ing Hamming distance as a metric. The problem is not 
rendered easy by the mere fact that the global minimum 
is so pronounced (although theoretical analysis inspired 
an efficient classical algorithm for the p = 2 Hopfield 
network, briefly described later on). It does imply, how¬ 
ever, that the ground state wavefunction does not jump 
chaotically: Every subsequent tunneling involves shorter 
distances, with 0{VN) spin flips, and achieves progres¬ 
sively better approximation to the true global minimum. 
Absent such a trend, annealing would be most difficult 
toward the end of the algorithm, when T ^ 1/A, and the 
minimum gap would exhibit less favorable exponential 
scalin^^. 

In what follows, the model and its solution are de¬ 
scribed in greater detail. First, finite-size scaling of the 
‘easy’ QCP bottleneck is linked to the thermodynam¬ 
ics of the phase transition. The next part goes beyond 
thermodynamics, considering small corrections to the ex¬ 
tensive contribution to the free energy. The entire low- 
energy spectrum, which depends on disorder realisation, 
is mapped onto that of a simple quantum mechanical 
particle in a random potential. Finally, extreme value 
statistics is applied to investigate the properties of that 
random potential near its global minimum by mapping 
it to a Langevin process. This yields the distribution of 
hard bottlenecks in a universal regime (T ^ 1). 

Quantum Hopfield network. Consider a model 
with rank-p matrix of interactions and no longitudinal 
field {hi = 0)P^ 

= ( 5 ) 

ii=i 

(cf. rk Jik = A for SK model), where are taken to be 
i.i.d. random variables of unit variance. The thermody¬ 
namics of this quantum Hopfield model has been worked 
out in great detail by Nishimori and NonomursP^. In 
fact, that study motivated the development of QA® 

When p is small {Jik ^ l/A^)? if is appropriate to 
replace local longitudinal fields with their m ean valu es 
The identity (erf) = hily/T"^ + hj is 
used to close this system of equations. For a transverse 
field below the critical value, T < Tc = 1, there appears a 
non-trivial (m ^ 0) solution to the self-consistency equa¬ 
tion for the macroscopic order parameter, a vector with 
p components: 

1 ^ 

(6) 

^ i = l 
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Here, the disorder variables are also written using vector 
notation: 4^^^). 

This model is equivalent to the Curie-Weiss (quantum) 
ferromagnet, which has a continuous phase transition 
characterised by a set of mean-field critical exponents. 
Two of these are particularly useful in the analysis of the 
annealing complexity: the one for the singular compo¬ 
nent of the ground state energy oc \j\^) as well 

as the dynamical exponent for the gap {Ei — oc 7^), 
where 7 = T — Tc is the ‘distance’ from the critical point. 
These exponents are defined for the infinite system, yet 
fairly general heuristic analysis (see the Methods section) 
predicts finite-size scaling for the QCP bottleneck: 

AEc(xN~'^^, ATc(xN~^^. (7) 

Substituting values a = 2 and b = 1/2 for the problem 
at hand, one may estimate the gap at the critical point 
and the width of the critical region to be and 

respectively. 

Worse-than-any-polynomial complexity of quantum 
annealing might be expected for the first order phase 
transition, which exhibits no critical scaling (but see 
ref. [33] for an exception). Another possibility is for the 
dynamical exponent to diverge at the infinite randomness 
QCP: the finite-size gap scales as exp(—c\/]V) in a ran¬ 
dom Ising chairP^. For the Hopfield model, however, this 
scaling is polynomial, as the disorder is irrelevant at the 
critical point. More intriguing is the fact that these pes¬ 
simistic scenarios are not found in SK spin glass either: 
the model is characterised by the same set of critical ex¬ 
ponents, albeit with logarithmic correct ion^^^M^. These 
corrections to scaling increase the gap and, respectively, 
decrease the width of the critical region by a factor of 
log^/^ N. Thus, as long as Tann ^ Af, non-adiabatic tran¬ 
sitions at the critical point should be suppressed. This 
presents a conundrum as SK model is known to be an NP- 
hard problem; finding a polynomial-time (even in typical 
case) quantum algorithm would be a surprising develop¬ 
ment. The heuristic analysis is clearly insufficient, but 
‘digging’ deeper into a problem would require a more 
‘microscopic’ analysis. In the following, the problem is 
mapped to ordinary quantum mechanics to uncover its 
low-energy spectrum that explicitly depends on a partic¬ 
ular realisation of disorder, 

Exact low-energy spectrum. Mean-field theory can 
be derived in a more systematic manner via Hubbard- 
Stratonovich transformation. Finite-temperature parti¬ 
tion function Z{/3) = Tr(e“^^) can be written as a path 
integral over m(t), which now acquires a dependence on 
the imaginary time 0 ^ ^ /d, with periodic bound¬ 
ary conditions: = m(0). The value of the inte¬ 

gral is dominated by stationary paths corresponding to 
the minimum of an effective potential V(m). While the 
discussion above has been deliberately equivocal on the 
distribution of disorder variables, it is now instructive to 
contrast bimodal = Al) Gaussian choices. The 
shapes of the effective potential for both scenarios are 


depicted in Fig. 



FIG. 2. The shape of the disorder-averaged effective po¬ 
tential V(m) below the critical point (F = 0.5 < Fc) for the 
choice of bimodal (left) and Gaussian (right) distribution of 
disorder variables. Red color highlights the minimum of the 
potential. For the bimodal distribution, the 2p-fold degener¬ 
ate global minima are organised in pockets corresponding to 
encoded patterns. For Gaussian distribution, the degenerate 
minima form a continuum, connected by arbitrary rotations. 

The conventional bimodal choice defines the model of 
associative memory: In the limit F = 0 the ‘patterns’ can 
be perfectly recalled {si = when p is small. For 

the Gaussian choice, the global minimum corresponds to 
a mixtur^^ 

Si = ±sgn(y]a^c7)> (8) 

rendering memory useless. In the bimodal case, such 
‘spurious’ states only become stable once the number of 
patterns scales with the problem siz^^ p > O.OSA/". The 
BCp (bimodal) or 0{p) (Gaussian) symmetry of the ef¬ 
fective potential is only approximate, to leading order in 
N. The degeneracy of the ground state is 2 (due to global 
spin inversion symmetry) for almost all disorder realisa¬ 
tions when p ^ 3 or p ^ 2 in the bimodal and Gaussian 
scenarios respectively. (Note that that the p = 2 bimodal 
case possesses an additional symmetry, which makes the 
ground state 4-fold degenerate.) The system is in a sym¬ 
metric superposition of the degenerate global minima at 
the end of QA. It evolves entirely in the symmetric sub¬ 
space since the time-dependent Hamiltonian commutes 
with U = exp(7ri5‘^). Thus, it should be noted that 
small gaps between symmetric and antisymmetric wave- 
functions are irrelevant to QA and can be ignored. 

Disorder fluctuations ‘nudge’ QA toward the ‘correct’ 
pattern as it passes the critical point in the bimodal Hop- 
field model. No further bottlenecks are encountered; gaps 
for F < Fc as well as the ‘classical’ (F = 0) gap are 0(1). 
By contrast, the classical gap scales as 0(1/^") in the 
Gaussian Hopfield model, alerting to a ‘danger’ posed by 
the F = 0 ‘fixed point’. To find the low-energy spec¬ 
trum when F < Fc, note that the dominant contribution 
to the path integral is from paths where the magnitude 
of the ‘magnetisation’ vector remains approximately con¬ 
stant, close to its saddle-point value, while the angle is a 
slow function of time: m(t) « {~coT'd(tf) P = 2. 
Integrating out the amplitude degrees of freedom, the 
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partition function is rewritten as 


Z{(3)<x j[dm]e 0^ ^ ^ , 


( 9 ) 


which describes a quantum-mechanical particle of mass 
M = 0{N) moving on a ring, subjected to a random 
potential 


Vrm = - ^ + N (v^)^ , 

( 10 ) 

where = 6( l^^e last term, written in short¬ 

hand, adds a constant offset so that {V{'3))^ = 0. Notice 
that Vr{'i^) oc \/]V via central limit theorem, thereby rep¬ 
resenting a higher-order correction to the extensive part 
of the free energy. 

Since the partition function Z{f3) = Q-I^En encodes 
information about the spectrum, low-energy (Goldstone) 
excitations of the many-body problem are in one-to-one 
correspondence with the energy levels of a quantum me¬ 
chanical particle, up to a constant shift. The next step 
is to find the properties of this potential near a global 
minimum. These are relevant in a regime away from the 
critical point, T ^ Tc, where the universal behavior char¬ 
acterised by the appearance of ‘hard’ bottlenecks sets in. 

Evolution of the random potential. Scaling of the 
gap for T < Tc can be obtained via semiclassical analysis. 
Small level splitting due to tunneling between wells at the 
two degenerate global minima (this degeneracy is a conse¬ 
quence of the global Z2 symmetry: Vr{'&-\-7r) = Vr{'&)) is 
not relevant to QA. Higher degeneracies are statistically 
unlikely; quantisation rules predict 0(A/'“^/^) gaps be¬ 
tween energy levels within the symmetric subspace. But 
this refers to the typical gap, obtained for fixed T for al¬ 
most all realisations of disorder. As quantum annealing 
sweeps the transverse field for a fixed realisation of disor¬ 
der, Vr{'d) might undergo global bifurcation. This would 
result in a small tunneling gap for a specific value of T 
when the competing minima are in resonance. 

Such a scenario is impossible near the QCP. Coeffi¬ 
cients in the Fourier expansion of the random potential, 
cos 2/c'd + 6/c sin2/c'd), decrease as so 

that the first harmonic dominates for T ^ Tc- Semiclas¬ 
sical analysis confirms a 0(A/'“^/^) gap at the critical 
point (where the curvature of the effective potential van¬ 
ishes, leaving only the quartic part). As T decreases, the 
random potential becomes more rugged (see Fig.|^ left), 
smooth only on scales A'd ^ F, which makes global bifur¬ 
cations more likely. Furthermore, it exhibits properties 
that allow one to make important predictions without de¬ 
tailed calculations. Rescaling the potential in the vicinity 
of either global minimum Vr{'d*) = Rp , 

( 11 ) 

V"r - Vr* ^ - ^r*), 

describes the same model but for the rescaled F ^ ^F and 
a different realisation of disorder. This scale invariance is 


responsible for the logarithmic scaling of the number of 
tunneling bottlenecks as has been explained earlier in the 
text. However, it still remains necessary to demonstrate 
that the density of bottlenecks is non-zero^ which entails 
an examination of the properties of the random potential 
in the limit F = 0. 




FIG. 3. Geometric interpretation of the random potential, il¬ 
lustrated for a random instance with A = 20. Vector m('d) is 
defined by drawing a separating line at angle ^ and assigning 
spins with on either side of separating line values = +1 
and Si = — 1 respectively (left figure). For changing 'd vec¬ 
tor m(^) is incremented/decremented by (2/A)^- (top right). 
The increments form a closed path, approximated by a circle; 
fluctuations define a random potential (bottom right). 


The classical optimisation problem corresponds to 
maximising the magnitude of ^ necessary con¬ 

dition for a local minimum is that two sets of vectors, 
{Ciki = +1} can be separated by a 

line. As the angle of this line changes, fluctuations of 
the amplitude give rise to a random potential Vb('d) (see 
Fig.§. This suggests a linear-time algorithm for finding 
a global minimum: Sort all vectors by angle (this may 
introduce an extra log N factor due to sorting overhead) 
and exhaustively check all possible angles. Of course, QA 
algorithm is too generic to exploit the specific structure 
of the problem; moreover ad hoc efficient algorithms are 
unlikely to exist for more general spin glass problems. 

On short intervals, the rando m pro cess is described 
as an undamped Langevin process^^^ in the continuous 
{N oc) limit (hence the exponent of 3/2 in eqn. 
corresponding to its fractal dimension). Properly tak¬ 
ing into account the statistics of the extremal proper¬ 
ties, the process must be conditioned on the fact that 
Vb('d) > Vb('^o) = ^0* away from the global minimum. As 
described in Methods and the Supplement, such a con¬ 
ditioned process consists of two uncorrelated branches, 

Vq{^) — Vq oc y/N I Integrating equations 

d(lnx±)/dr = z/±(r), 

Aq/A -X- 
did/dr = ±x± , 

defines X+('d) and X-('^) parametrically, in terms of ran¬ 
dom processes that correspond to a Brownian mo¬ 

tion in a non-linear potential depicted in Fig. This 
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FIG. 4. Illustration of appearance of local minima. Left: 
Potential Vr{0) for a specific realisation of disorder as a func¬ 
tion of r. A perspective-projection 3D plot that zooms on a 
region near the global minimum of Vq{0) is shown in the bot¬ 
tom. Right: Top part of figure plots a specific realisation of a 
stochastic process z/(t) and a Lan gev in potential U.(y) (inset). 
The bottom part integrates eqn. ( p^ to obtain x('r) and ^(r) 
in the lower figure. Lower inset plots x as a function of'd. 
Regions with z/(t) < 0 are responsible for the appearance of 
local minima. 


potential is biased toward positive ‘velocities’ v so that 
Vq —V"o* ^ from eqn. (12). It will, however, 

experience arbitrary percentage drops due to subpaths 
with < 0 (albeit with decreasing probability). 

For small but finite F, the ‘separating line’ becomes 
blurred. The random potential adds a ‘quantum correc¬ 
tion’ (see Methods and the Supplement), Vp {'if) — Vo {if) = 


0(\/]vr3/2). For two minima to come into resonance, 
they cannot be more than Aif ^ F apart (i.e. 0{NT) 
spin flips). The tunneling exponent is given by under- 
the-barrier action A ^ VMAVAif^ where the effective 
mass M ^ N/T‘^ and AV ^ VTvr^/^, leading to eqn. 
Numerical results for the universal distribution of tunnel¬ 
ing jumps and the tunneling exponents are exhibited in 
Fig.|5l 


DISCUSSION 

Poor scalability of classical annealing of spin-glass 
models had been linked to the phenomenon of temper¬ 
ature chaoP^. Interestin gly, it s existence in mean-field 
glasses had been debatecP^^®, alt hough it is uncontro- 
versial in finite-dimensional model^^^^. Similarly, fail¬ 
ures of quantum annealing might be attributed to trans¬ 
verse field chaos. The phenomenon described in this 
manuscript represents a much stronger finding. The mere 
fact that ground states at F and F + AF will have vanish¬ 
ingly small overlap as A ^ cxo is not inconsistent with the 
continuously evolving ground state and poses no ‘threat’ 
to QA. By contrast, ‘hard’ bottlenecks correspond to iso¬ 
lated discontinuities that persist as AF ^ 0. 

To dwell upon the generality of these results, first note 
that scaling of the tunneling exponent will depend on 


A0/r 
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FIG. 5. Numerical results. The L-shaped figure (composed 
of panels I-III) plots the distribution of the tunneling jumps 
and the exponent. The lower left quadrant (panel I) is a 
scatterplot of A'd/T (proportional to the number of spin flips 
normalised by TN) against | log AF^tunnl/ {TN)^^^ (the prefac¬ 
tor c) for all tunneling events. Histograms of these quantities 
alone (projections of the scatterplot) are in the bottom right 
(panel II) and the top left (panel III) quadrants. The top 
right quadrant (panel IV) plots a histogram of distribution of 
Fn/Fn+i for the successive tunneling events. 


the universality class of the model. The SK model, for 
instance, exhibits different scaling of barrier heights, be¬ 
lieved to be oc (see e.g. ref. |44]). In the model 

studied, the decrease in the number of spins involved in 
tunnelings offsets the divergence of the effective mass in 
the classical limit F ^ 0. As the height of the barriers 
also decreases, the tunneling gaps widen toward the end 
of the algorithm. One might expect qualitatively similar 
behavior in realistic spin glasses. 

The logarithmic scaling of the number of bottlenecks 
is due to self-similar properties of the free energy land¬ 
scape in the interval [Tmin^rc]- The lower cutoff should 
correspond to the smallest energy scale in the classical 
limit, which for the SK model is also a negative power of 
A, namely l/v^A- This is related to the linear vanishing 
of the density of distribution of effective fields as /z —> 0 
at zero temperatur^^ (since P{h)dh ^ 1/A). The 
picture is less clear for constraint satisfaction problems, 
where the energies are constrained to be non-negative in¬ 
tegers; i.e. the classical gap is 0(1). These energy levels 
have exponential degeneracy, which is lifted by the trans¬ 
verse field. A value sufficient to make the spectrum ef¬ 
fectively quasi-continuous might serve as an appropriate 
lower cutoff Fmin in problems of this type. Lack of ‘hard’ 
bottlenecks in the Hopfield model with the bimodal dis¬ 
tribution of disorder and p = 0(1) could be attributed 
to the fact that the number of valleys is finite, which is 
not representative of a true spin glass. 

An interesting observation is that since ‘hard’ bottle¬ 
necks correspond to Landau-Zener crossings, annealing 
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times need not be exponentially small. The probability 
that QA fails to follow the ground state every single time 
in n repeated experiments is 

(e-5(TO)”, (13) 

which exactly matches the probability of failure for the 
annealing rate that is n times slower. Using shorter an¬ 
nealing cycles with many repetitions can minimize the ef¬ 
fects of decoherence. It suffices that non-adiabatic transi¬ 
tions are suppressed at the critical point only, ^ < 1/N. 

Even with polynomial annealing rates, coherent evolu¬ 
tion would require much better isolation from the envi¬ 
ronment than what is currently feasible. The only com¬ 
mercial implementation of QA (D-Wave) must contend 
with a fairly strong coupling to a thermal bath. On the 
positive side, it accommodates faster annealing cycles, 
acting as a ‘safety valve’ to dissipate any heat generated 
during the non-adiabatic process. At the same time, it 
all but ensures that the system is always in thermal equi¬ 
librium with the environment. 

The theory presented here breaks down when T < T so 
that quantum bottlenecks described here may not be a 
limiting factor if ln(rc/T) < l/a. Main source of errors 
will be from exponentially many thermally occupied ex¬ 
cited states. If the annealing profile were adjusted so that 
the energy spacing increased beyond T towards the end, 
this would effectively implement classical annealing. An 
intricate relationship between temperature, problem size, 
and the properties of the spin glass model might deter¬ 
mine which mechanism (quantum or classical annealing) 
will be dominant. 

Yet another tradeoff in the design of D-Wave chip is a 
quasi-2D topology of interactions Jik due to fabrication 
constraints, which incurs significant performance penalty 
when mean-field models are ‘embedded’ into a ‘Chimera’ 
grapl]^. So-called ‘native’ problems corresponding to 
uniformly random instances on this quasi-2D lattice have 
been argued to be poor candidates for QA due to the 
lack of a finite-temperature classical phase transitioiP^, 
at the same time, a quantum phase transition at Tc > 0 
is expected in 2D glasseJ^. 

Whereas l®^-order phase transition immediately im¬ 
plies exponential complexity, even for small sizes, prob¬ 
lems having a continuous phase transition may remain 
tractable up to a threshold, Nc^ beyond which tunnel¬ 
ing bottlenecks become dominant {alnNc ^ 1). This 
‘tractability threshold’ serves as a silver lining fot this 
otherwise negative result. Moreover, the picture of ‘hard’ 
bottlenecks may be equally applicable to classical anneal¬ 
ing. While in some crafted examples classical anneal¬ 
ing is at a unique disadvantage due to -order phase 
transitiorP^, for most interesting problems both classical 
and quantum transitions are 2^^-order. In such a sce¬ 
nario, the density of bottlenecks becomes a tie-breaker 
for evaluating relative performance. Whether quantum 
annealing can be advantageous in terms of this metric 
and determining which models will benefit is a practi¬ 
cally important question for follow-up work. 


There remains a question whether the failure mecha¬ 
nism described here can be somehow circumvented. Such 
a feat is feasible, e.g. for a disordered Ising chain, where 
an exponentially small gap develops at the critical point, 
which is a manifestation of Griffiths singularities^. Mod¬ 
ification of QA that requires controlling the transverse 
field for each spin individually can suppress these singu¬ 
larities and restore a polynomial gap. 

Based on comparison with another exactly solvable 
model, it seems that frustration — in addition to dis¬ 
order — is essential for the appearance of small gaps in 
the spin glass phase proper. The seemingly random pro¬ 
file of the energy landscape for finite T heralds difficul¬ 
ties in avoiding these bottlenecks in generic spin glasses. 
Although any prospects of exponential speedup should 
be met with skepticism, heuristics inspired by spin glass 
theory revolutionalised branch-and-bound algorithms for 
constraint satisfaction problem^^. One can remain hope¬ 
ful that theoretical advances can similarly aid quantum 
optimisation. 


METHODS 

Finite-size scaling of the gap at QCP is best understood using an 
example of a finite-dimensional system. In thermodynamic limit, 
both correlation length and characteristic time diverge near the 
phase transition as 

oc (r-Tc)-’', Tc oc (14) 

In a finite system this divergence is smoothed out as soon as the 
correlation length becomes comparable to lattice size (^c ~ 

The minimum gap (the reciprocal of r) and the width of the critical 
region should scale as N~^l^ and respectively. In this 

paper, the product zv has been labeled as exponent b. Singular 
behavior of normalised ground state energy (free energy) is related 
to the specific heat exponent (a = 2 — a). Dimensionality d can be 
eliminated with the aid of hyperscaling relation 2 — a = (d + z)^ 
to yield eqn. 0 of the main text. Independent estimates of the 
specific heat exponent can be obtained from the exponents for the 
order parameter and the susceptibility (2 — a = 2/3 + 7 ). 

Finite-temperature partition function can be written as a sum 
over a set of paths [sdOl with 0 ^ t ^ /3, where Si(t) alternates 
between the values ±1. Hubbard-Stratonovich transformation can 
be used to rewrite it as a path integral 

e5w7(Ei€iSi(t))''dt+i:;Kr[si(t)] 

(15) 

~ j [dm(l)]e-f fo inZi, 

The ‘kinetic’ term JCr[s(t)] in the first equation penalises kinks, 
representing F-dependent ferromagnetic couplings between Trotter 
slices. As the interaction term is decoupled, the problem reduces 
to that of evaluating the single-site partition function Zi for a spin 
subjected to a magnetic field with the transverse component F and 
the ‘time-dependent’ longitudinal component hi(t) = ^^ m(t) . To 
leading order in A, the saddle-point of the path integral ( |15| > is a 
solution of mean-field equations. This becomes a degenerate man¬ 
ifold for Gaussian disorder distribution; to determine higher-order 
contribution all paths such that \m(t)\ = mr are considered. Eval¬ 
uating Zi is best performed in the adiabatic basis that diagonalises 
the 2-level Hamiltonian Hi = —hi{t)d'^ 

Zi = TrTe^o = Tr Te“ 


(16) 



Here Ei(t) is diagonal with eigenvalues Its fluctu¬ 

ations around the mean give rise to the random potential 
Non-adiabatic terms due to rotation of the basis {U = are 

treated using 2nd order perturbation theory, giving rise to a kinetic 
term oc (d-d/dt)^. Note a simple form of the perturbative term in 
eqn. ( |S7| ) owing to the fact that Vi(t) commute for all t. 

In the limit N —)■ oo, the random potential defined in 

eqn. in the main text converges to a Gaussian process that can 
be specifled completely by its covariance matrix (Vr('*^)hr'('*^0)- 
This can be ‘diagonalised’, alternatively expressing the random 
potential as a linear combination of white-noise processes {Cn('d)}. 
One representation, as a convolution series (/r^^^Cn) ('d) with 

kernels 

oc ^r2 + ( 2^2 ^ (17) 

relies on orthogonal properties of associated Laguerre polynomials. 
The choice o; = 1 ensures that only n = 0 term survives in the limit 
r = 0. The random potential should satisfy a stochastic equation 


dVo 




(IS) 


As a side note observe that a similar equation is obtained by taking 
a continuous limit in the identities that follow from elementary 
geometry (see Fig.j^in the main text, 


Vk+l-Vk _ Vk-Vk-1 ^ 
cfe+1 Cfc 

+4) + 4y214/Arcos cos ( 19 ) 

For finite F, the form of the potential is modified as follows: 
It is convolved with a smoothing kernel of width A-d ~ F, which 
raises the global minimum by 0(\/iVF^/^). Additional random 
contributions (from n ^ 1) have similar scaling. 

In the vicinity of the global minimum, the statistics of the clas¬ 
sical potential is fundamentally different. The Returning force’ in 
eqn. ( |18| ) can be neglected; additionally the process should be con¬ 
ditioned on the fact that it stays above its value at -d = 0 (no gen¬ 
erality is lost by choosing the global minimum as the origin). This 
problem has been a subject of a considerable body of worl®I^, 
although important aspects ought to be revisited. Here, I present 
a particularly simple self-contained derivation. 


A pair (x? '^) — where x cx: Vq — is interpreted as the ‘coor¬ 
dinate’ and V = dx/d'd as the ‘velocity’ (-d being ‘time’) — forms 
a Markov process. The probability distribution p('d; x, '^) satisfles, 
for -d > 0, 

^ = + while p(i?;+0, d) = 0 (V-u > 0) (20) 

serves as a boundary condition for the absorbing boundary. This 
probability is ‘renormalised’ to condition on the fact that it sur¬ 
vives until some © ^ -d. It becomes a conserved quantity, but the 
diffusion equation adds a drift, —d(\og Pq)/ dv^ repelling from the 
boundary. The ‘survival’ probability P© in the limit 0 —)► oo is, 
up to ‘time-reversal’, the universal asymptotic solution of ( |20| ), re¬ 
duced to ordinary differential equation (ODE) using scaling ansatz: 


p('d, Xj '^) oc 


x2a/3^^ (^/xl/3) 


( 21 ) 


This exploits a fact that fractal dimensions of ‘velocity’ and ‘coor¬ 
dinate’ are [i;] = and [x] = [ 7 ?]^/^. The asymptotics is dom¬ 

inated by solutions with the smallest possible exponent, a = 1/4 
out of the inflnite set of eigenvalues for the ODE. This matches a 
known value obtained with a different methocS^. 

The next step performs a change of variables, introducing ‘di¬ 
mensionless’ velocity u = and ‘logarithmic’ coordinate 

= Inx- With the ‘time’ variable redeflned via dr = d'd/x^'^^, 
Markov process is described by a tuple ( 7 ?, fi, u). Marginalising out 
11 and 'd in the equation for p(t; t?, fi^ u) produces Fokker-Planck 
equation for a stochastic motion of particle in a potential 


U{v) = - In ((i//6^/®) Ai(i/2/6^/®) + Ai'{v'^. (22) 


Given a particular realisation of z^(t), the full process Hq i/, 7?) is 
determined deterministically, by integration (see eqn. ( |12| ) in the 
main text). The construction of a realisation of a random process 
is performed independently for 7 ? > 0 and 7 ? < 0. 

Numerical simulations rescale this random potential instead of 

_ AT _ 3 AT 

evolving the transverse held: 7 ? 1 —)■ e r -d and x ^ 2 r y. The 
process is extended to larger values of r as needed (details of the 
process for small r, where they fall below the numerical precision, 
are ‘forgotten’). A fairly large range of r is required to gather the 
sufficient statistics. 

The code used for simulations and the raw numerical results are 
available upon request. 
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SUPPLEMENTARY INFORMATION 


Mathematical Preliminaries. This paper makes extensive 
use of Meijer G-function defined via intergal 


Qm,n (^ ni, . . . , Un; ^n+l, • • • , \ _ 

V bl, . . . ,bm',bm+l, .. ■ ,bg J 
± f nr=l + «) r(l -at-s) 

^^iJcUl=r.+inak + s)UUr.+ini-bk-s) ’ ^ > 


where the contour C is chosen appropriately. Bold T(x) is used 
to denote E-function to avoid confusion with the transverse field 
variable. G-function can be used to express many elementary and 
special functions; those used here include: 

• (1 + x) “ = 0 )’ 

• Lii \x)e ^ 

• U{a,b,x) = r(a)r(a-6+i)^i’2(^l o,i-b)’ 

where L^\x) are associated Laguerre polynomials and U{a,b,x) 
represents confiuent hypergeometric (Kummer) function. Mul¬ 
tiplication by a finite power x^ merely shifts all parameters 
ai,..., an, ^1, • • •, by constant a. Most integrals encountered 
can be recast as those of Mellin type: 





p+r,g+s 


a, c\ 

b, d)’ 


(S2) 


where vectors are concatenated as follows: a, c consists of ai,..., Un 
and Cl,..., C£ followed by Un+i ,. .. ,ap and C£_^i,. .., Cr (similarly 
for b and d). This ensures that zeros and poles in the Mellin trans¬ 
form of G-functions are grouped together. 

Mean-Field Analysis. Self-consistency equation for the ‘mag¬ 
netisation’ vector defined in eqn. § of the main text can be written 
as follows: 


m = 


i 


^ i \/r2 + (4im)2 ’ 


(S3) 


where the individual spins polarise along the magnetic field with 
the components E and hi = in the transverse and longitudi¬ 
nal directions respectively. Replacing sum over spins with disorder 
average, one obtains non-trivial solution to self-consistency equa¬ 
tions for E < 1. For the bimodal distribution, the solutions are 
(±\/l — E^, 0,..., 0), (0, ±\/l — E^,... , 0), etc. Other spurious so¬ 
lutions appear for smaller E, but they do not become stable for 
finite p. For Gaussian distribution, the magnetisation vector can 
have arbitrary direction while the magnitude mr satisfies 


1 = 


/ 








where the integral on the right hand side had been converted 
to Mellin type, | ^o^)^o’i(^| 7 ) ^ substituting t = 

^/2 and x = r^/( 2mp). Alternatively, it can be expressed in 
terms of modified Bessel functions. 

Low-energy spectrum of Hopfield model can be obtained by ex¬ 
amining the partition function at finite temperature. It can be 
written as a sum over paths [si(t)] alternating between values 
Si(t) = ±1 with periodic boundary conditions Si{(3) = Si(0): 


e5 So Kr[si(t)], 


(S5) 


The ‘kinetic’ term /C[s(t)] = (# of kinks) x ^ IncothTAt, where At 
is the discretisation chosen (limit At —)► 0 will be taken eventually). 
This term completely suppresses the dynamics in the limit E = 0. 


With Jif^ = (l/A') two-body interactions are de¬ 

coupled via Hubbard-Stratonovich transformation: Using the iden¬ 
tity oc f dme“^"^ the partition func¬ 

tion is recast as a path integral, 

y'[dm(t)]e-f(gg) 

where 

— 'y ^ qJq hi(t)si(t)dt-\-)Cp[s(t)] _ (S7) 

[s{t)] 


is a single-site partition function for a spin in a ‘time-dependent’ 
longitudinal field hi{t) = ^^m(t). It has been rewritten with the aid 
of time-ordering operator T. The single-site Hamiltonian Hi{t) = 
-hi{t)cr^ -Vcr^. 

The dominant contribution to the path integral is given by 
stationary paths, m(t) = m. To leading order In this case 
Z rkj J where the effective potential is 


V(m) = - (^r2 + (Cm)2)^ = - V^|m|C/(-|;0; ^). 

(S8) 

The saddle-point value mr = |m| that minimises the potential for 
r = 0 is given by a solution of eqn. ( |S4| ). 

Gonsidering higher-order corrections, the dominant contribution 
comes from paths where the magnitude of magnetisation fiuctuates 
around the mean-field value and the angle is a slow function of 
time: m(t) mr( local field hi{t) is slow- 

varying, it is convenient to evaluate the partition function ( |S7| > 
in the adiabatic basis that diagonalises the instantaneous 2-level 
Hamiltonian: Hi(t) = where 

Eiit) = — and Vi{t) = - arccot ^ . (S9) 

Using this factorisation, the time-ordered product 

can be rewritten in the limit At —)■ 0 as 

Zi =TrTe~ (SIO) 

The non-adiabatic term can be simply written as i-^V{t) since 
all V{t) oc Q-y are commuting (should be —U^^U more gener¬ 
ally) and acts as a perturbation. The lower-energy level, E^(t) = 

—+ h? (t), makes a dominant contribution to the partition 
function in the limit /3 —)■ oo. Including second-order correction 
from the perturbation theory, \dV^f /dt\‘^/(E^ — E^)^ one obtains: 


\nZi 



r^{dhi/dtf \ 

8(r2 + /l2)5/2^ 


(Sll) 


As the sum over all sites is performed, the first term in the integrand 
gives rise to a finite-size correction to the effective potential, which 
depends on angle 7? and a particular realisation of disorder 


Vp(,?) = - ^r2 + sm2(i? - di) + 

(S12) 

Vector-valued has been parametrised here as Off¬ 

setting average in eqn. ( |S12| ) above has been made with respect to 
Gaussian distribution of the projection of ^ onto m (its value is 
indep ende nt of the direction). The second term in the integrand of 
eqn. pil l contributes 

Jo 2 Y4(r2+$2^2gi„2(^_0.))5/2 

^ ^ ^ 

M 
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The ‘effective mass’ M can be replaced by its average value 


M = 


Nmr 

4V2r2 


1 

2 ’ 


- 1 ; 


r2 


)■ 


(S14) 


which approaches N/ in the limit F —)■ 0. The dynamics 

becomes more classical as the effective mass diverges. Path integral 
( |S6| ) is then rewritten as 

z ^ e-^'3V(mr) I [di»(t)]e“ (S15) 

The low-lying energy levels of the many-body problem are thus 
equivalent to those of a quantum mechanical particle on a ring. 

The approximations made above break down for any finite N 
when F is sufficiently small. To avoid this, and to be able to make 
more q uanti tative predictions about the behavior of random po¬ 
tential ( |S12| , it is necessary to study the continuous limit N ^ oo 
directly. 

Continuous Random Potential. Observe that the cen¬ 
tral limit theorem can be applied to the entire two-dimensional 
(in F and 7?) process Vr(^)- This Gaussian process can be 
decorrelated by writing it as an infinite function series = 

En m Cnm/nm(r5 involving independent Gaussian variables 
Cnm- Basis functions are determined from the covariance 
. Since this correlation depends only on the angle 
difference -d — -d', it is natural to use /nm(F, 7?) = /nm(F) 

Although the best truncated approximation is obtained by us¬ 
ing Karhunen-Loeve basis, there is no requirement that /nm(F) be 
orthogonal. In fact, since covariance involves averaging over dis¬ 
tribution V{^) = it is quite convenient to use a set of 

polynomials that are orthogonal with said weight. Using a set of 
white noise processes 


{CnWCn'i'^')) = ^nn' 7^,7?' G [-7r/2; 7r/2], (S16) 


This matches the correct value obtained by replacing sum over sites 
by disorder averages and determines the normalisation constant: 


A-n, — 


2^7rr(7T, + q; + 1) 


(S21) 


Notice that the value of a can be chosen freely: as mentioned 
above, the decomposition is not unique. The most convenient 
choice, a = 1, leads to the vanishing of all convolution kernels 
with 72 ^ 1 in the limit F = 0. Introducing the rescaled potential 
via Vb('d) = (2\/N /7r)x('d), it is straightforward to verify that 


xW= ||sin('d - 6>)| Co(6>)d6> (S22) 

satisfies the following stochastic equation: 

S+X = CoW. (S23) 

The solution is determined uniquely by enforcing periodic boundary 
conditions: x(-7t/ 2) = x(7^/2) and x'(-7r/2) = x'(7^/2). 

Now considering a finite F > 0, it becomes convenient to re¬ 
express 72 = 0 term in the functional series as a convolution with 
x('d) instead. Thi s can be a ccom plished by substituting the left 
hand side of ( |S23| > into eqn. i |sT7| i and performing integration by 
parts. Other terms can also be more conveniently expressed by 
rewriting white-noise processes in terms of continuous processes: 
realisations of periodic brownian motion 

^+Cn=CnW. (S24) 

dv 

Here Cn — I- ^2Cn('d)d7? is added to ensure periodicity: 

77n(-7r/2) = 77n(7r/2). 

Now, the random potential is rewritten as 


in lieu of discrete random variables, the series expansion becomes 
°° r^/2 , . 

VrW = VNj2 4"\^-e)<:„{e)de, (si7) 

where a factor \/iV had been introduced for convenience. The con¬ 
volution kernels should be chosen so as to match the covariance. 
One suitable choice is 


(S18) 

Since associated Laguerre polynomials form a complete orthogonal 
(with respect to weight F^e“^ with t = /2) set of basis function, 

it follows that 


oo 


2'^r(72 + a + I)An / 2 2 C2 • 2 q 

-^^-a/F^ + 772^^^ sm^ 7?. 


(S19) 

Representing Ur('*^) using ansatz of eqn. ( |SI7| ) above one can write 
the expression for the covariance matrix. Using identities ( |SI8| > 
and ( |S19| >, it is transformed as follows: 


oo „ 7L 
n=0 "'-f 

oo „ 7L 

= NAn 

n=0 "'-f 

J ^r'2 +7712,^2 sin2(^/ _ 


2°‘r{n + a + l)NAl 

nl 



F^ + 7T2p^2 sin^ {'d — 6) X 


^r'2 + m2,^2 sm2(,?' (g2o) 


. (S25) 


The last term is angle-independent random variable that absorbs 
{fn}- Being a constant offset, it cannot affect the physics; and, 
thus, will be neglected. Expressions for kernels fr('*^) and are 
obtained by performing integration by parts in variable 7?: 



F2 (f 2+^2^^2 

2(F2 + ^2 gjj ^2 ^^^/2 






sinT^cosT? 

2VnTT Jo 


(526) 

(527) 


where the rescaled transverse field F = F /TUp has been intro¬ 
duced for convenience. Using substitutions t = ^/2 and x = 
r^/(2sin^ 7?), the integrals above can be converted to Mellin form 
and evaluated in terms of special functions: 


frW 





3 ^2 \ 

2’°’ 2sh?tfy 


+ t'^u 


5 

2 ’ ’ 2 sin^ 7? ) 


COST? ^2,2/ 

\/n + 1 n! ^’®\2sin^i? Ii2;0 ) 


(528) 

(529) 


Tunneling gaps are associated with global bifurcations of Up ('d) 
which occur mostly in the limit F <C I. Moreover, it is this limit 
that is governed by universal scaling laws and thus will be investi¬ 
gated in greater detail. The kernels above act on scales 7^ ^ F . For 
small F it is permissible to drop the second term in eqn. ( |S28| > and 
to replace sin 'd 'd and cos 7 ? ps I everywhere. It will be seen that 
any bifurcations also take place for 7? ~ F, hence the periodicity of 
the potential Up(7?) becomes irrelevant. 

Extremal Statistics. The distribution of bottlenecks in the 
limit F <C I is related to the properties of the classical potential 
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xi'^) in a proximity of its global minimum It will be governed 
by different statistics, derived below, by virtue of this extremality 

condition. _ 

On short scales, the linear term on the left hand side of eqn. ( |S23| ) 
can be neglected in comparison with the large stochastic term 
Co ('*?)• It will be seen self-consistently that the relevant range is 
I'd — t^qI = 0(r). The stochastic equation describes a free parti¬ 
cle in ID subjected to random force. Specifying the ‘coordinate’ 
X together with the ‘velocity’ v = dy/d'd for some d determines 
their distribution either in the ‘future’ (d' > d) or in the ‘past’ 
(d' < d), where d plays the role of time. This can be described by 
the following PDE: 


= 0, 


(S30) 


dp dp 1 d^p 

d'd dx 2 dv"^ 

where p(d; x? '^) is the probability density. It will be convenient to 
perform a shift so that a global minimum coincides with the origin: 
dQ = 0 and x(0) = 0- Another condition is necessary to ensure 
that x('^) > 0 lor d / 0. Considering only d > 0 for concreteness, 
one may introduce an absorbing boundary via boundary condition 


lim^p(d; x,'T') = 0 for 'L’ > 0. (S31) 


Notice that the probability need not vanish as x —^ +0 for negative 
velocities 'i; < 0. The probability that a random path does not hit 
this boundary decays as follows: 

^ J p(d;x,^)dxd'u = J i;p(d;+0, i;) d-u. (S32) 

Since the random process is conditioned on the fact that it starts at 
a global minimum, it is convenient to ‘renormalise’ the probability 
so that it becomes a conserved quantity once again. Consider a 
probability that x stays positive at least until some distant horizon 

B, 


dr = dd/x^'^^j if is now a triple (d, x, that becomes new Markov 
process. This modificatio n mer ely adds + f ('r’/x)^ fo 

the left hand side of eqn. ( |S35| . The term oc q reflects renormalisa¬ 
tion of probability density as a result of non-linear transformation. 

Furthermore, it is convenient to introduce new ‘dimensionless’ 
velocity u = vjx^^^ and the ‘logarithmic’ coordinate p = \nx- The 
new probability density, appropriately renormalised, reads 

q(T] d, p, v) = e^^/3^(T; d, e^, e^/^). (S38) 

This is substituted into eqn. ( |S35| > which has been modifled as de¬ 
scribed in the previous paragraph. Replacing djdx {d/dp — 

^vd/dv) and djdv ^ e~^/^dldu, the following PDE is obtained: 


dr d'd dp, 



1 d'^q 
2^ 


= 0, 


(S39) 


where the non-linear potential is written, using eqn. ([S37|, as 


W(0 ~ Y -lnp.(-0- 


(S40) 


When p(r; d, /r, v') is marginalised over d, p, the standard Fokker- 
Planck equation is obtained. It describes a stochastic process 


dz/ 

dr 


-U'(v) sgnT + C(r), 


(S41) 


where C(^) is the white-noise random force. This stochastic equa¬ 
tion can be integrated both forward (r > 0) and backwards 
(r < 0) starting from z/(0) drawn from the equilibrium distribu¬ 
tion p{u) oc 

To And the shape of U(y) and p{y )^ ob serve that p*(y 
satisfles a stationary version of eqn. ( |S30| >, which is reduced to 
ordinary differential equation. Writing 2 ; = 


gf(d; X?'^) oc p(d; Xj'^) [ ^(B; X, T|d; x,'^) dX dT, (S33) 

^©>0 




with an appropriate normalisation factor. The second factor, 
Pq{'&; XiP) represents the survival probability. It is writte n as an 
integral of the Green’s function associated with eqn. ( |S30| ). Time- 
reversal symmetry implies 


P(©; X, T|d; X, v) = P{'^; X, -'^|B; X, -T). (S34) 

As a result P©(d; x, p) satisfles a PDE obtained from eqn. ( |S30| ) by 
changing the sign of ^ and With this in mind, the equation 

for renormalised probability becomes 


dq dq d / 1 dPQ \ I d‘^q 

d'd ^X ^ \ Po dv ^) 2 dv‘^ 


(S35) 


In comparison with eqn. ( |S30| ), it adds a drift oc d{\ogP<Q)/dv in 
addition to the diffusion of ‘velocity’ v. 

In the limit of large 0, the ‘survival’ probability P©(x,'r’,'d) 
corresponds — u p to time-reversal ('L’ —)■ —'L’) — to the asymptotic 
solution of ( |S30| , 

(S36) 


Probability density always converges to this form as 'd ^ 00 , 
with dependence on the initial conditions only via the prefactor 
A. Here p*(X;' ^) ca n be thought of as a ‘stationary’ solution: 
it satisfles eqn. ( |S30| > ( with d/d'd term dropped) and boundary 
condition ( |S31| ). Eqn. ( |S30| is invariant with respect to rescal¬ 
ing (-d, X, W ^ (M, pi ‘^v). For this reason, the asymptotic 

solution should be sought in the form 

P* (x, v) = P* • (S37) 


It is possible to simplify equations even further by introduc¬ 
ing new ‘dimensionless’ variables. Deflning new ‘time’ variable via 


z'lp" + ^'0' + = 0, (S42) 


where 'ip{z) = a/ p{z). Disallowing exponentially increasing solu¬ 
tions, 'ip{z) is sought in the form e^^f±{z) for z>_0 and 2 ; < 0 
respectively. This substitution transforms eqn. ( |S42| > to hyperge¬ 
ometric form, giving the two branches as f-\-{z) oc U §5 22 ;) 

and f-(z) oc P (| + ^; ^;—2z). Matching logarithmic deriva¬ 
tives 


^>(+0) r(-^) 


and 


V>(-o) r(| + f) 


(S43) 


[where the common factor is A = 3^/^r(|)/r(|)] requires that 


a = 


1 

4 


3 

+ -n, 


(S44) 


as verifled using of Euler’s identity r(a:)r(l — x) = tt/ sin(7ra:). 
Of this infinite set, the smallest positive value a = 1/4 has to be 
chosen. The analytical expression that deflnes the potential U{v') 
can be simplifled as follows: 

'ip{iy) oc = X where x = v j v^. (S45) 

Since eqn. ( |S39| > is first order in 7 ?, these are completely deter - 
ministic, given a particular realisation of stochastic process ( |S41| : 

Y = 

^=X^/^(r), (S47) 


with boundary conditions x(~c)o) = 'd(— 00 ) = 0. This procedure 
deflnes the process x('^) parametrically. So far, a positive branch 
(-d > 0) had been considered. Negative branch (-d < 0) is obtained 
similarly; positive and negative branches are uncorrelated. 









